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Abstract 

Implementations of solid state quantum optics provide us with devices where qubits are placed at fixed 
positions in photonic or plasmonic one dimensional waveguides. We show that solely by controlling the po- 
sition of the qubits and with the help of a coherent driving, collective spontaneous decay may be engineered 
to yield an entangled mesoscopic steady-state. Our scheme relies on the realization of pure superradiant 
Dicke models by a destructive interference that cancels dipole-dipole interactions in one-dimension. 
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The study of atoms coupled to the electromagnetic (e.m) field confined in cavities or optical 
waveguides has played a central role in the fields of Quantum Optics and Atomic Physics. In recent 
years the basics of that physical system has been realized with artificially designed atoms in solid- 
state setups. We may include here quantum dots and nitrogen vacancy (NV) centers determinis- 
tically coupled to photonic cavities [1-3], and plasmonic [4, 5] or photonic [6-10] waveguides, 
as well as circuit QED setups where superconducting qubits are coupled to microwave cavities 
[11, 12]. Even though the physics of atomic and solid-state quantum optical systems is similar, 
the latter shows a crucial advantage: on a solid substrate, emitters may be placed permanently at 
fixed positions at separations of the order of relevant wavelengths [13]. An important application 
of those systems is quantum information processing in solid devices [14], where artificial atoms 
acting as qubits are placed within the e.m. field confined in a microcavity. Typically, the realiza- 
tion of those ideas requires unitary qubit- field evolutions and a natural candidate is to consider 
collective couplings to a single mode in a cavity. An alternative pathway is to tailor the interaction 
with the environment to induce quantum correlations between qubits with dissipation [15]. This 
approach has been proven to be advantageous to generate entanglement between ensembles of 
atoms [16]. In this direction one-dimensional guided modes have been recently pointed out as an 
useful tool to create two-qubit entanglement [17-19] and many-qubit entanglement in cascaded 
quantum networks [20]. 

In this Letter, we show that by placing a set of qubits in a one-dimensional waveguide (see 
Fig. 1) the continuum of e.m. field modes induces a controllable dissipative coupling between 
the qubits. The possibility of deterministically position the artificial atoms or qubits allows us to 
engineer the paradigm for quantum optical collective effects, i.e. the Dicke model of superradiance 
[21] in its pure form. The observation of the latter in optical systems is hindered due to dephasing 
caused by dipole-dipole interactions [22, 23]. In our scheme those interactions can be switched 
off by an appropriate choice of the inter-qubit distance. Adding a classical drive to the pure Dicke 
model we obtain a dissipative system with a phase diagram of steady states showing mesoscopic 
spin squeezing and entanglement. This model has been theoretically investigated in the past [24- 
26], but experimental realizations are scarce. Finally we upgrade our scheme to a set of N 4-level 
emitters [27, 28] and show that a judicious choice of couplings to the waveguide and dispersion 
relations may lead to a variety of many-body dissipative models which show entangled steady- 
states. 

We start by modeling N 2-level systems (2LS), {\g) n , \e)n}n=i...N> placed at positions x n and 
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Figure 1. (Color Online) (a): Experimental scheme of the system: ensemble of equally- spaced qubits 
placed in the vicinity of a one-dimensional waveguide, (b): Two-level system configuration with resonant 
excitation. Panel (c): Four-level system configuration with two-additional lasers. 

coupled to a one dimensional field (Fig 1 (a)) with photon annihilation operators a q , described by 
the Hamiltonian H — H^+Hi. The free term is Hq — H q b + Afield > with (h = 1) 

#qb = -y L G ni Afield = £ ® q a\a q , (1) 

where C0o is the qubit energy (Fig 1 (b)) and co q is the field dispersion relation. We define Pauli 
matrices G% = \e) n (e\ — \g) n (g\, <*t = \ e )n(g\> G n — \s)n{e\- The photon polarization is neglected 
to focus on the most relevant physics of our work. We consider a dipolar coupling of the form 

#i = £(a„£(x„)+H.c.), (2) 

n 

with E{x) = Y, q gq{ a q eiqx + a\e~ iqx ), and g q is a dipolar qubit-field coupling. Define p, the re- 
duced density matrix for the qubits. In the weak coupling limit, the evolution of p can be described 
by a markovian master equation of the form dp jdt — J£f (p) [29], with the superoperator 

^(P) = K"P^m " PO-+CJ-) +H.C. (3) 
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A detailed derivation follows the description in dimensions higher than one presented in previous 
works [23] (see Section A Sup. Info). Special care must be paid to the counter rotating terms in 
(3), which have to be included to get the following result for the collective decay rates 

/ — _Jq(coo)\x n -x m \ / 4 x 

We define F = y(coo), with the function y(o)) = g 2 ,,D((o)/K, where q(co) is the resonant 
wavevector at 0), and we have defined the e.m. density of states, D(co) = (2k / L)\dq(co) / dco\, 
with L the quantization length. A crucial observation for this work, is that couplings J n ^ m do not 
decay with the distance [30], a situation that is singular of one dimensional waveguides. In free 
space on the contrary collective couplings decay like 1/r or 1/r 3 , depending on the relative dipole 
orientation [23]. 

Homogeneous couplings [17, 18, 31] / n?m — T/2 can be obtained from Eq. (4) by the choice 
x n = ftAo, with Ao = 2k /qo, and n £ Z. This condition cancels dipole-dipole interactions and we 
get the pure Dicke superradiant decay described by 

JS? D (P) = \ (S-pS + -S + S-p) +H.c, (5) 

with S~ = Y,n a n> s+ =Hn G n- We also define S a = o„/2 (a = x,y,z) 9 and the basis {\J,M)} 
of eigenstates of S 2 , S z . Assuming an initial state like |*Fo) = ®n\e)n — \N /2,N /2) the system 
evolves within the sector J = N/2. We note that Dicke superradiant decay is achieved in one 
dimension without the restriction that the whole atomic ensemble is confined within a region of 
length Ao. 

In this work we focus on the qubit steady-state, p s , which for a given Liouvillian fulfills 
J£?(p s ) = 0. To achieve some controllability on p s , we add a pump term which physically can 
be implemented by the interaction of qubits with a resonant field with Rabi frequency £2, 

&>*(p)=&>(p)-i%[S x ,p]- (6) 

Competition between the collective decay and the pumping leads to a non-equilibrium phase tran- 
sition in the steady-state of the model at a critical pumping rate £l c = AT/2 [24], manifested in 
a kink in the population inversion observable (S z ), see Fig. 2(a). Let us first give a brief de- 
scription of the two limiting cases: (i) Coherent steady state regime, £1 <C NT/ 2. Since ^b,p 
can be obtained from J^b by the substitution S~ —> S~ + /£2/(2r), one can easily show that 
p s = pF c )(¥ c | + <^ 2 (^), where |¥ c ) = e i T s *\N/2,-N/2) is a spin coherent state, (ii) Mixed 
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Figure 2. (Color online) Numerical results for the coherently pumped Dicke model (6). Different color 
represent different number of qubits from N = 2 (blue) to N = 200 (red) as shown in the legend, (a) 
Population inversion (S z )/N. (b) Purity & (dashed) and spin- squeezing parameter l/^ 2 (solid). 

state phase, £1 3> NT/ 2. Here we get an infinite temperature state. To show it, it is conve- 
nient to write J??d in the interaction picture with respect to £lS x /2. This accounts to replace 
S~ — >> S x + (1/2) (cos (^Sy + sw(t)S z ). Averaging over time, leads to 

J?Vv^(s x pS x -S 2 xP + l - £ (SapSa-Slp^j+H-C. (7) 

which has the infinite temperature state p s = 1 as steady-state. For calculations in the intermediate 
regime we use the full solution in the |/,M) basis. 

To quantify many-particle entanglement we use spin- squeezing £ as a figure of merit, 

2 N(AS y f 
g (S x ) 2 +(S z r W 

The latter is both an entanglement witness and it is also linked to applications in quantum metrol- 
ogy [32-34]. Values £ < 1 are spin-squeezed states which can be proved to be entangled [33]. 



Note that the above mentioned (i) and (ii) phases lead to £ = 1 and £ = 0, respectively. Another 
theoretical tool to be used is the purity, defined by £P(p) — Tr(p 2 ). Both magnitudes are plotted 
for different number of qubits N in Fig. 2. We find a range of pumping fields (£1 < £l c ) that in- 
duce a pure entangled p s . This result may lead to the generation of many-particle correlations in 
mesoscopic samples of artificial atoms. 

We upgrade now our 2LS into a 4-level system (4LS) configuration, see Fig. 1 (c), and 
show that this is a more advantageous situation. Our scheme can be realized in the solid-state 
context [27, 28] and describes a variety of possible configurations in which a set of low-level 
states are coupled to excited states by lasers with different polarizations. Two ground states 
(|#)nj \ e )n) are coupled to high energy states (\g f ) n , \e') n ). The qubit part of the free Hamiltonian 
becomes now — Ln( fi Vl£%(£ , | + ®A e, )n{e'\ + (O g \g) n {g\). Two weak non-resonant fields 
with amplitudes ^l,i(2) an d frequencies G)l,i(2)> induce transitions described by a Hamiltonian 
term// L = Y, n {{^\/2)\e') n {e\e- i ^ t + (£l U2 /2)\g f ) n (g\e- i(OL ^ + H.c). After an adiabatic elim- 
ination of the excited states (see Sec. II of Supp. Mat. for details) we get an effective qubit-field 
interaction H\{t) — H^ q (t) +Hf(t), written in the interaction picture with respect to Hq. The first 
term reads 

Hi\t) ^KEMiDle^+Dne-^'), (9) 

n 

where K 2 = (£2l,i/2Ai) 2 — (0,^2/ 2A 2 ) 2 (Ai^) = GV(g') ~~ ^,1(2)) * s a normalization constant. 
D n = ua~ +v<7+ is a jump operator resulting from the cross radiative decay, with u = k* -1 £2l,i/2Ai, 
v = k* -1 £2l,2/2A2, fulfilling u 2 — v 2 = 1. The latter condition will be useful in the discussion be- 
low, and allows us characterization by a single parameter r, such that u = cosh(r), v = sinh(r). 
After eliminating the e.m. field degrees of freedom we arrive to the Liouvillian 

^sq(p) = £/* q m {D-pD+-pD+D-) +H.c. (10) 

with f^ m = r sq ^^(^)l Xw_Xm l and T sq = K 2 y((O a ). The second term in the effective qubit-field inter- 
action describes longitudinal decay processes, 

= KY^E(x n ,t)(ua^e iC0lt +va z e- i0>2t +U.c.), (11) 

n 

and leads to 

^z(p) = L<m(^p^-P<^)+H.c 5 (12) 
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with/^ m = K 2 (y(coi)u 2 e iq ^ G)l ^ Xn ~ Xm \ + y^a^v 2 e iq ^\ Xn ~ Xm \) . J£ z induces a dephasing mechanism 
that competes with the spontaneous coherence build up induced by j£f q . The relative importance 
of those contributions depends on the photon density of states at frequencies 0) a and (D\ 2. We 
consider two limiting cases: (i) Small photonic bandwidth.- This is the most favorable configura- 
tion. We assume that the density of states in the waveguide is peaked around 0) a , with a bandwidth 
Aco <C \(0\ — G) a |, |g>2 — &>a| such that 7(6)1,2) ~ and therefore J* « 0. For example, this can be 
the case of one dimensional waveguides consisting of coupled cavities forming a one-dimensional 
photonic crystal. Defining q(co a ) = 2n/?i a and choosing x n = nX a we arrive to a spin- squeezing 
version of the Dicke superradiant model, 

^sq, D (p) = (D-pD + -D + D-p+H.c.) , (13) 

where we have introduced collective squeezed spin operators D + /~ = E w AT^ + . In Fig. 3 (a) 
we present a calculation of the spin- squeezing in the steady- state as function of the squeezing 
parameter, r. Remarkably, we observe an enhancement of the maximum value of entanglement of 
several orders of magnitude compared to the case of an ensemble of 2LS's. (ii) Large photonic 
bandwidth.- In the opposite limit we consider a broadband waveguide [17, 31] (|cDl,i — <*>l,2| <^ 
Aco) such that the density of states at the frequencies considered here is comparable y(co\) ~ y(coi) . 
In experiments with optical transitions, for example with quantum dots in optical or plasmonic 
waveguides, condition C0\, CO2, co a 3> coo is found, since transition energies are in the eV and meV 
for high energy and low energy transitions, respectively [27, 28]. Thus, we can safely assume 
q(co a ) & q(C0i) ~ q(coi) = 2n /X a , and consider that quantum dots can be placed at the same 
relative optical path with respect to all frequencies. To give a more quantitative argument for 
this approximation we define the group velocity of the modes of the waveguide v g (co) = \d q co q \, 
and consider the limit \co\ — a) a |v g (a) a ), \a>i — 0) a |v g (cOb) <C g(ft) a ), which corresponds to small 
wavevector differences. In the case of constant v g and optical transitions, this condition leads to 
differences of 10 -3 in ^(o) a ), q{(0\^)- We neglect for the moment those differences, which may 
lead to inhomogeneous broadening effects that are discussed later in this work. Thus, he condition 
x n — nX a leads to a collective dephasing term of the form 

^z,D(p) = y (S z pS z - S z p + H.c.) , (14) 

where we have introduced the rate T z — J^ n . In the large photonic bandwidth limit we get thus two 
competing terms Jzf = J2? S q,D +^fz,D- Collective dephasing increases with the squeezing parameter 
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r, as depicted in solid black in Fig. 3 (b). The competition between dephasing and squeezing 
mechanisms determines an optimal r to generate maximal entanglement. The latter can be higher 
than the one generated by «£b, p considered above. The purity of the system is also affected by the 
dephasing term, however one can still find a region that combine high purity and high values of 
entanglement as shown in Fig. 3 (b). 

Finally, we discuss the feasibility of our ideas focusing on the following points: (i) One- 
dimensional waveguides.- We require long propagation length and efficient coupling to the qubits. 
Significant coupling to the guided modes of more than 85 — 89% have already been reported for 
photonic [9, 10] and plasmonic [4, 5] waveguides. Theoretical predictions of even higher efficien- 
cies has also been pointed out [31, 35] at expense of increasing the losses and reducing the prop- 
agation length. Limitations in the propagation length will limit in practical terms the number of 
qubits one can entangle. Besides, the precise location of the qubits is also needed, which is possi- 
ble for solid-state emitters using, i.e, lithographic methods which nowadays have precision bigger 
than 50nm [13]. (i) Lambda-transitions in solid-state qubits.- We assume a degree of manipulation 
of electronic levels similar to the one achieved in atomic physics, specially the 4LS scheme. Appli- 
cations in quantum information processing [14, 36] typically require controlling optical transitions 
for spin-pumping and initialization. Recent experimental results [27, 28] show level schemes in 
quantum dots similar to those required in our work. (Hi) Markovian approximation.- We require 
FN <C C0o in the 2LS scheme and FN <C G)li,2? in the 4LS case, such that the cooperative de- 
cay rate is much smaller than the transition frequencies, the latter determine the photonic bath 
memory time [29]. This condition is well satisfied in the case of optical transitions of quan- 
tum dots, (iv) Independent decay channels on each transition frequency.- This is required for 
the 4LS's scheme in the large bandwith limit, to get Eqs. (10, 12), and it is justified as long as 
FN <C | 0>l 5 i — C0 a | ? | <*>L,2 — I • Since differences in transition energies are of the order of meV , this 
condition imposes a restriction on the achievable rates for entanglement generation in our scheme. 
(v) Homogeneous couplings.- Our scheme neglects inhomogeneous couplings which would take 
the steady state out of the collective basis |/,M). This may be a severe restriction since inhomo- 
geneous broadening is a limitation of quantum optical solid-state devices. Assuming a dispersion 
of energies induced by a term = AcjjSj, we expect a dephasing time l/t^oc var(Afi)/), such 
that entanglement is destroyed for times t > T^. We have checked this result with exact numerical 
calculations for small size systems. Although inhomogeneous broadening in solid-state setups is 
still of the order of meV for quantum dots [13] and jieV for NV centers [37], the feasibility of our 
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proposal will benefit from current experimental efforts in the field. 

In conclusion, we have proved that one-dimensional plasmonic [4, 5] or photonic [6-10] waveg- 
uides are excellent candidates to correlate a large number of qubits by collective radiative decay. 
By exploiting the site-controlability of solid state emitters [13] we get an exciting playground 
to study non-equilibrium dissipative dynamics. We have presented several schemes with two or 
four level systems which are ideally suited to achieve mesoscopic entanglement in controllable 
solid-state devices currently under investigation for quantum information processing. A thorough 
analysis on the experimental challenges has been done. Our ideas can be easily translated to circuit 
QED by generating of sidebands in the qubit-field coupling [38]. 

During completion of this work we became aware of a related theoretical preprint on the prepa- 
ration of spin-squeezed states in atomic ensembles [39]. Acknowledgments.- We acknowledge 
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Figure 3. (Color online) Same as Fig. 2 for the 4LS scheme, (a) Limit of Jf nn = in the entanglement 
witness l/£ 2 as a function of the parameter r for different numbers of qubits. (b) Entanglement witness 
(l/£ 2 ) and purity (^) as a function of the squeezing parameter for different number of qubits. 
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